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Abstract. This paper presents a review of the phenomenon of interlayer ex- 
change coupling in magnetic multilayers. The emphasis is put on a pedagogical 
presentation of the mechanism of the phenomenon, which has been successfully 
explained in terms of a spin-dependent quantum confinement effect. The theo- 
retical predictions are discussed in connection with corresponding experimental 
investigations. 
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1. Introduction 

Since the first observation by Griinberg et al. (lj of antiferromagnetic coupling 
of Fe films separated by a Cr spacer, the interlayer exchange interaction between 
ferromagnetic layers separated by a non-magnetic spacer has been a subject of intense 
research in the last few years. The decisive stimulus came from the discovery, by 
Parkin et al. [Q, of oscillations of the interlayer exchange coupling in Fe/Cr/Fe and 
Co/Ru/Co multilayers, as a function of spacer thickness. Furthermore, Parkin || 
showed that this spectacular phenomenon occurs with almost any metal as a spacer 
material. 

This review will be restricted to the case of interlayer coupling across non- 
magnetic metallic spacer layers. This excludes the cases of non-metallic spacer 

S pl Js, || i jjl Pll, of antiferromagnetic spacer such as Cr or Mn [|Io|, |l2| , |l3| |l4| , 
|16|, |l7|, |J, jig |(| |l| and of rare-earth multilayers ||, ||[ |i[j25). This choice is 



motivated by the fact that the physical mechanism of the coupling in these cases is 
quite different from the one to be discussed here. 

The magnetic coupling energy per unit area can usually be expressed as 

E(9) = Jcos9, (1) 

where 9 is the angle between the magnetizations of the two magnetic layers, and J is 
called the interlayer coupling coupling constant. Higher order terms in an expansion 
in powers of cos 9 are also observed; such terms, which give rise to non-collinear 
alignment of the magnetizations, are believed to be of non intrinsic origin and to be 
related to defects such as roughness ]2(| ^7j. Such effects will not be considered here. 
In addition, due to space limitation, the important question of the role of alloy disorder 



and interdiffusion 28, E9l [M BlJ will not be addressed. 
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The purpose of this paper is to present as simply as possible the mechanism 
of interlayer exchange coupling in terms of quantum interferences due to electron 
confinement in the spacer layer. The understanding of this mechanism relies on ideas 
due to various authors (3^, gjj, |34|, |35|, p6fl . The presentation give here is based on 
the one of Ref. |56|. The emphasis will be on physical concepts and pedagogical 
clarity rather than on mathematical rigor. It is organized as follows: in the next 
Section, an elementary discussion of quantum confinement is given; in Section |3], it is 
then shown how spin-dependent confinement in the spacer layer gives rise to interlayer 
exchange coupling; Section [| is devoted to the limit of large spacer thicknesses, for 
which particularly simple results are obained; Sections |E| and ^ treat the variation of 
interlayer exchange coupling with magnetic layer thickness and non-magnetic overlayer 
thickness, respectively; finally, in Section [?] the strength of the interlayer exchange 
coupling is discussed in comparison with experimental data. 

The point of view adopted here reflects the author's subjective views on the 
topic. In particular, due to space limitation, the important literature devoted to ab 
initio calculation will not be discussed in detail here. The interested reader can find 
complementary information in the various review papers on this subject which have 
been published recently [§§ H g ||| [|| || g§ . 

2. Elementary discussion of quantum confinement 

For the sake of clarity, we shall first consider an extremely simplified model, namely 
the one-dimensional quantum well, which nevertheless contains the essential physics 
involved in the problem. Then, we shall progressively refine the model in order to 
make it more realistic. 

The model consists in a one-dimensional quantum well representing the spacer 
layer (of potential V = and width D), sandwiched between two "barriers" A and B 
of respective widths La and Lb, and respective potentials Va and Vb- Note that we 
use the term "barrier" in a general sense, i.e., Va and Vb are not necessarily positive. 
Furthermore, the barrier widths, La and Lb, can be finite or infinite, without any 
restriction. 

2.1. Change of the density of states due to quantum interferences 

Let us consider an electron of wavevector k + (with k + > 0) propagating towards 
the right in the spacer layer; as this electrons arrives on barrier B, it is partially 
reflected to the left, with a (complex) anplitude rs = |rB|e 1( ^ B . The reflected wave 
of wavevector k~ is in turn reflected on barrier A with an amplitude ta = I^Ale 1 ^" 4 , 
an so onjj] The module |rA(B) | of the reflection coefficient expresses the magnitude of 
the reflected wave, whereas the argument 4>a(b) represents the phase shift due to the 
reflection (note that the latter is not absolutely determined and depends on the choice 
of the coodinate origin). 

The interferences between the waves due to the multiple reflections on the barriers 
induce a modification of the density of states in the spacer layer, for the electronic 
state under consideration. The phase shift resulting from a complete round trip in the 
spacer is 

Ac/) = qD + <$> A + 4>b , (2) 

| Of course, for the one-dimensional model, one has k~ = —k + ; however, this property will generally 
not hold for three-dimensional systems to be studied below. 
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with 

q = k+ -k~ . (3) 

If the interferences are constructive, i.e., if 

Acp = 2mr (4) 

with n an integer, one has an increase of the density of states; conversely, if the 
interferences are destructive, i.e., if 

Acj> = (2n + 1)tt (5) 

one has a reduction of the density of states. Thus, in a first approximation, we expect 
the modification of the density of states in the spacer, An(e), to vary with D like 

An(e) re cos (qD + (t> A + 4>b) ■ (6) 

Furthermore, we expect that this effect will proportional to the amplitude of the 
reflections on barriers A and B, i.e., to |r>tPfl|; finally, Aro(e) must be proportional to 
the width D of the spacer and to the density of states per unit energy and unit width, 

7r de 

which includes a factor of 2 for spin degenaracy. We can also include the effect of 
higher order interferences due to n round trips in the spacer; the phase shift A(f> is 
then multiplied by n and Ir^re is replaced by It'a^sI™- Gathering all the terms, we 
get, 

An(e) « — ^ V \r A r B \ n cos n (qD + 4> A + fa) 
7r de 

n— 1 

iDf £ ±(r A r B r e^A 

n—1 / 

4 r ^ qD D ) (8) 
de 1 -r A r B e lqD J 

As will appear clearly below, it is more convenient to consider the integrated density 
of states 

N{e)= ^ n(e')de'. (9) 

J — oo 

The modification A7V(e) of the integrated density of states due to electron confinement 
is 




AN(e) = UmJ2 {l ^-e niqD 



11 

n—l 

2 ^ - '- jqD\ 



= Im In (l — r A r B c lq ) (10) 

7T 

A simple graphical interpretation of the above expression can be obtained by noting 
that Im ln(z) = Arg (z), for z complex; thus, AN(e) is given by the argument, in the 
complex plane, of a point located at an angle A(f> = qD + <f) A + <f> B on a circle of radius 
Ir^r^l centred in 1. This graphical construction is shown 
The variation of A7V(e) as a function of D is shown 
the confinement strength r^r^l. For weak confinement (a), A7V(e) varies with D 
in sinusoidal manner. As one increases the confinement strength (b), the oscillations 



in |IL 
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for various values of 



Theory of interlayer exchange interactions 



4 



Im(z) 



i 














/V B I \ 








Re(z) 


Arg(z) \ 


1 J 





Figure 1. Graphical interpretation of equation (|lo|). 



are distorded, due to higher order interferences. Finally, for full confinement (c), 
AN(e) exhibits some jumps that correspond to the appearance of bound states. We 
note however, that the period A of the oscillations of AN(e) does not depend on the 
confinement strength, but only on the wavevector q = k + — k~ , namely, A = 2ir/q. 

So far, we have implicitely restricted ourselves to positive energy states. Negative 
energy states (i.e., of imaginary wavevector) are forbidden in absence of the barriers 
A and B, because their amplitude diverges either on the right hand side or on the left 
hand side, so that they cannot be normalized. This matter of fact no longer holds in 
the presence of the barriers if Va (or Vb, or both Va and Vb) is negative: the negative 
energy states, i.e., varying exponentially in the spacer, can be connected to allowed 
states of A or B. In order to treat these states consistently, we simply have to extend 
the concept of reflection coefficient to to states of imaginary wavevector, which is 
straightforward. One can check that, with this generalization, ( |Io| ) acounts properly 
for the contribution of the evanescent states. Physically, this can be interpretated as 
a coupling of A and B by tunnel effect |36f.. 



2.2. Energy associated with the quantum interferences in the spacer 

Let us now study the modification of the energy of the system due to the quantum 
interferences. In order to conserve the total number of electrons, it is convenient 



to work within the grand-canonical ensemble, 
grand-potential, which is given by 



and to consider the thermodynamic 



4> 



In 



1 + exp 



ef — £ 



n(e) de 
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Figure 2. Variation of AN(e) as a function of D, for various values of the 
confinement strength: (a) Ir^r^l = 0.1, (b) |rA?\B| = 0.8, (c) r^rg = 1 (full 
confinement). Note the different scales along the ordinate axis. 



N(e)f(e)de. (11) 



At T = 0, this reduces to 

(e - ef) n(e) Ae 

-co 

= - f F N(e)de. (12) 

J — oo 

The energy AE associated with the interferences is the contribution to $ 
corresponding to AN(e), 

9 r+oc 

AE=-lm In (l - r A r B e iqD ) de. (13) 
2.3. Three-dimensional layered system 

The generalization of the above discussion to the more realistic case of a three- 
dimensional layered system is immediate. Since the system is invariant by translation 
parallely to the plane, the in-plane wavevector k|| is a good quantum number. 
Thus, for a given k||, one has an effective one-dimensional problem analogous to 
the one discussed above. The resulting effect of quantum intereferences is obtained 
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by summing on k|| over the two-dimensional Brillouin zone. The modification of 
integrated density of states per unit area is 

1 T f ,9. , /. jq±D\ 



AN(e) = - — Im / d 2 k|. In (1 - r A r B e^ v ) , (14) 

and the interference energy per unit area is 

1 f f +co 

AE=^lmj d 2 k|, J /(e) In (l - r A r B e^ D ) de . (15) 

2.4- Quantum size effect in an overlay er 

The case of a thin over layer deposited on a substrate is of considerable interest. In this 
case, one of the barriers (say, A) consists of the vacuum, and barrier B is constituted 
by the substrate itself. The potential of the vacuum barrier is V vac = ep + W, where 
W is the the work function; thus it is perfectly reflecting for occupied states, i.e., 
\fvac \ = 1- On the other hand, the reflection on the substrate (or coefficent r su b) may 
be total or partial, depending on the band matching for the state under consideration. 

The spectral density of the occupied states in the overlayer can be investigated 
experimentally by photoemission spectroscopy; in addition, by using inverse 
photocmission, one can study the unoccupied states. If furthermore these techniques 
are used in the "angle-resolved" mode, they give information on the spectral density 
locally in the k|| plane. 

For a given thickness of the overlayer, the photoemission spectra (either direct or 
inverse) exhibit some maxima and minima corresponding, respectively, to the energies 
for which the interferences are constructive and destructive. When the confinement 
is total, narrow peaks can be observed, which correspond to the quantized confined 
states in the overlayer, as was pointed out by Loly and Pendry flUjfl , 

Quantum size effects due to electron confinement in the photoemission spectra of 
ove rlay ers have been observed in various non-magnetic systems Q, ^5, 46[ 47], ^8|, [|9|, 
[sol , plj HI). In particular, the systems Au(lll) /Ag/vacuum and Cu(lll)/Ag/vacuum 
offer excellent examples of this phenomenon J49|, pl| . 

2.5. Paramagnetic overlayer on a ferromagnetic substrate: Spin-polarized quantum 
size effect 

So far our discussion concerned exclusively non-magnetic systems. Qualitatively new 
behavior can be expected when some of the layers are ferromagnetic. A case of 
particular interest is the one of a paramagnetic overlayer on a ferromagnetic substrate. 

In the interior of the overlayer, the potential is independent of the spin; therefore 
the propagation of electrons is described by a wavevector k± which is spin-independent. 
The reflection coefficient on the vacuum barrier, r va c, is also spin-independent. 
However, the ferromagnetic substrate constitutes a spin-dependent potential barrier; 
thus, the substrate reflection coefficients for electrons with a spin parallel to the 
majority and minority spin directions of the substrate, respectively r,! ub and 
It is convenient to define the spin average 



sub CL±AI - 1 '"sub 



sub 1 sub (16) 



and the spin asymmetry 



T _ 1 
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In this case, the electron confinement in the overlayer gives rise to a spin-dependent 
modulation of the spectral density versus overlayer thickness; the period of the 
modulation is the same for both spins, whereas the amplitude and phase are expected 
to be spin-dependent. 

The quantum size effects in paramagnetic overlayers on a ferromagnetic substrate 
have been investigated by several groups |5^, [54|, |55|, |5(| |57], p£| |5{| [H], ^ 
p4] , |65[ The systems studied most are Cu overlayers on a Co (001) substrate 

and Ag overlayers on a Fe(001) substrate. Ortega and Himpsel |54|, |35| observed a 
quantum size effect in the normal-emission photoelectron spectra of copper overlayer 
on fee cobalt (001) substrate. They observed peaks due to quantum size effects both 
in the photoemission and in the inverse photoemission spectra an oscillation of the 
photoemission intensity. These quantum size effects manifest themselves also by an 
oscillatory behavior of the photoemission intensity at the Fermi level; as the observed 
oscillation period (5.9 atomic layers) is close to the long period of interlayer exchange 
coupling oscillations in Co/Cu(001)/Co, they suggested that the two phenomena 
should be related to each other; they also claimed that the observed oscillations in 
photoemission are spin dependent and due mostly to minority electrons. A direct 
confirmation of this conjecture has been given independently by Garrison et al. [fi7f 
and by Carbone et al. |58|| by means of spin-polarized photoemission. They found 
that both the intensity and the spin-polarization exhibit an oscillatory behavior with 
the same period (5-6 atomic layers), but have opposit phases, which indicates that 
the quantum size effect does indeed take place predominantly in the minority-spin 
band as proposed by Ortega and Himpsel |34], pq] . Recently, Klasges et al. Q 
and Kawakami et al. ]66[ ] have observed spin-polarized quantum size effects in a 
copper overlayer on cobalt (001) for a non-zero in-plane wavevector corresponding to 
the short period oscillation of interlayer exchange coupling in Co/Cu(001)/Co; they 
observed short-period oscillations of the photoemission intensity in good agrement 
with the short-period oscillations of interlayer coupling. This observation provides a 
further confirmation of the relation between quantum size effects in photoemission 
and oscillation of interlayer exchange coupling. 

Photoemission studies of quantum size effects have also been performed in other 
kinds of systems such as ferromagnetic overlayer on a non-magnetic substrate, or 
systems comprising more layers [^37], |(|| [nj, [nj . 

Photoemission spectroscopy undoubtly constitutes a method of choice for 
investigating quantum size effects in metallic overlayers: this is due to its unique 
features, which allow selectivity in energy, in-plane wavevector, and spin. 

Besides photemission, spin-polarized quantum size effects in paramagnetic 
overlayers on a ferromagnetic substrate are also responsible for oscillatory behavior 
(versus overlayer thickness) of spin-polarized secondary electron emission [73fl , 
linear (74[ || || f||j7§ and non-linear [^0[ [H] magneto-optical Kerr effect, and 
magnetic anisotropy p2 , |83| . However, these effects usually involve a summation over 
all electronic states, so that the quantitative analysis of the quantum size effects may 
be fairly complicated. 



3. Interlayer exchange coupling due to quantum interferences 

Let us now consider the case of a paramagnetic layer sandwiched between two 
ferromagnetic barriers A and B. Now, the reflection coefficients on both sides of 
the paramagnetic spacer layer are spin dependent. A priori the angle between the 



Theory of interlayer exchange interactions 



8 



magnetizations of the two ferromagnetic can take any value; however, for the sake of 
simplicity, we shall restrict ourselves here to the ferromagnetic (F) configuration (ie., 
9 = 0) and the antiferromagnetic (AF) one (i.e., 9 = 7r). 

For the ferromagnetic configuration, the energy change per unit area due to 
quantum interference is easily obtained from flq), i.e., 



1 

47T 3 



Im 



d 2 k 



-foe 



In 1 



' A 1 S e 



In 1 



r A r B e 



de 



(18) 



In this equation, the first and the second term correspond respectively to majority- and 

minority-spin electrons. The antiferromagnetic conguration is obtained by reversing 

t i 

the magnetization of B, i.e., by interchanging r B and r B ; thus the corresponding 
energy per unit area is 



AE AF = — Im 

47T 3 



In 1 



d^k 



+ OC 



r A r B e 



In 1 



r l p i<?j_-D 
r A r B e 



de . (19) 



Thus, the interlayer exchange coupling energy is 



E A i 



1 

47T 3 



Im / d 2 k 



which can be simplified as 
E F - E AF 



x In 



1 



T 



A B 



r l r l p iq±D 
' A r B t} 



r A r B e 



)(" 



r A r B e 



de 



Im / d^k 



/(e) Ar A Ar B e iq±D de 



(20) 



(21) 



in the limit of weak confinement. The above expression for the IEC has a rather 
transparent physical interpretation. First, as the integrations on k|| over the first two- 
dimensional Brillouin zone and on the energy up to the Fermi level show, the IEC is a 
sum of contributions from all occupied electronic states. The contribution of a given 
electronic state, of energy e and in-plane wavevector k|| , consists of the product of three 
factors: the two factors Ata and Ars express the spin-asymmetry of the confinement 
due to the magnetic layers A and B, respectively, while the exponential factor e lq±D 
describes the propagation through the spacer and is responsible for the interference 
(or quantum size) effect. Thus, this approach establishes an explicit and direct link 
between oscillatory IEC and quantum size effects such as observed in photocmission. 



4. Asymptotic behavior for large spacer thicknesses 

In the limit of large spacer thickness D, the exponential factor oscillates rapidly with e 
and k|| , which leads to substantial cancellation of the contributions to the IEC due to 
the different electronic states. However, because the integration over energy is abruptly 
stopped at sp, states located at the Fermi level give predominant contributions. Thus 
the integral on e may be calculated by fixing all other factors to their value at e^, and 
by expanding q± = fc^ — fcj around e^, i.e., 

qx nq XF + 2^f^, (22) 
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U ±_F V ±F V ±F 

The integration (see Ref. |36) for details) yields 

E F -E AF = ^ Im | d 2 k„ Mi Ar A Ar B e^' D 

x F{2nk B T D/hv+ F ), (24) 

where 

= (25) 

In the above equations, q± F is a vector spanning the complex Fermi surface; the 
velocity v±_p is a combination of the group velocities at the points (k||,£fjj) and 
(k||, k^ F ) of the Fermi surface. 

Next, the integration on ky is performed by noting, that, for large spacer thickness 
-D, the only significant contributions arise from the neighboring of critical vectors kjf 
for which q±p is stationary. Around such vectors, q± F may be expanded as 

a {k x -Kf (ky-k*) 2 

q±F = q ±F — — , (26) 

where the crossed terms have been canceled by a proper choice of the x and y axes; 
k" and Ky are combinations of the curvature radii of the Fermi surface at (k|, k±) 
and (kf,A;- Q ). 

The integral is calculated by using the stationary phase approximation | f36"| , and 
one obtains 

E F - E AF = Im £ gg|ArSAr B e^ 

a 

x F{2nk B TD/hv c l), (27) 
where g", v", Ar^, Ar B correspond to the critical vector k", and 

KqS K) 1/2 ( K «) 1/2 ; (28) 

in the above equation, one takes the square root with an argument between and tt. 

This analysis shows that in fine, the only remaining contributions in the limit 
of large spacer thickness D arise from the neighborhood of states having in-plane 
wavevectors k? such that the spanning vector of the Fermi surface q± F = k^ F — k^ F is 
stationary with respect to ky for ky =k^, and the corresponding contribution oscillates 
with a wavevector equal to q^_ F - This selection rule was first derived in the context of 
the RKKY model it is illustrated in Fig. | There may be several such stationary 
spanning vectors and, hence, several oscillatory components; they are labelled by the 
index a. 

The above selection rule allows to predict the oscillation period(s) of the interlayer 
exchange coupling versus spacer thickness just by inspecting the bulk Fermi surface of 
the spacer material. In view of an experimental test of these predictions, noble metal 
spacer layers appear to be the best suited candidates; there are several reasons for this 
choice: 

• Fermi surfaces of noble metals are known very accurately from de Haas-van 
Alphen and cyclotron resonance experiments p4| ; 
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Figure 3. Sketch showing the wavevector giving the oscillation period of 
oscillatory interlayer exchange coupling, for the case of a non-spherical Fermi 
surface. 



• since only the sp band intersect the Fermi level, the Fermi surface is rather simple, 
and does not depart very much from a free-electron Fermi sphere; 

• samples of very good quality with noble metals as a spacer layer could be prepared. 

Fig. U shows a cross-section of the Fermi surface of Cu, indicating the stationary 
spanning vectors for the (001), (111), and (110) crystalline orientations [[33|; the Fermi 
surfaces of Ag and Au are qualitatively similar. For the (111) orientation, a single 
(long) period is predicted; for the (001) orientation, both a long period and a short 
period are predicted; for the (110) orientation, four different periods are predicted 
(only one stationary spanning vector is seen in figure |], the three others being located 
in other cross-sections of the Fermi surface) . These theoretical predictions have been 
confirmed successfully by numerous experimental observations. In particular, the 
coexistence of a long and a short period for the (001 ) orientation has been confirmed 
for Cu @ H, H, || |3, Ag H, and Au g |c], |j; and the experimental periods 
have been found to be in excellent agreement with the theoretical ones. A comparison 
of the theoretically predicted oscillation periods and the experimentally observed ones 
is given m Table @. 

In a further attempt to test the theoretical predictions for the periods of oscillatory 
coupling, several groups |)6[ |97| , |98| have undertaken to modify in a controlled manner 
the size of the Fermi surface (and hence, the period of the coupling) by alloying the 
spacer noble metal (Cu) with a metal of lower valence (Ni); in both cases, the change in 
oscillation period due to alloying has been found in good agreement with the expected 
change in the Fermi surface. 



Theory of interlayer exchange interactions 



11 




Figure 4. Cross section of the Fermi surface of Cu along the (110) plane passing 
through the origin. The solid dots indicate the reciprocal alttice vectors. The 
dashed lines indicate the boundary of the first brilouin zone. The solid arrows, 
respectiveley horizontal, oblique, and vertical, indicate the vectors giving the 
oscillation period(s), respectively for the (001), (111), and (110) orientations. 



5. Effect of magnetic layer thickness 

As already mentioned, the influence of the IEC on the ferromagnetic layer thickness is 
contained in the reflection coefficients Ar^ and Ar^. If the ferromagnetic layers are 
of finite thickness, reflections usually may take place at the two interfaces bounding 
the ferromagnetic layers, giving rise to interferences p9| , and hence, to oscillations 
of the IEC versus ferromagnetic layers thickness. A more detailled discussion of this 
effect is given in Refs. [B6L 99|. This behavior was first predicted from calculations 



based upon a free-electron model [ 100 1 . The amplitude of the oscillations of the IEC 
versus ferromagnetic layers thickness is generally much smaller than the oscillations 
versus spacer thickness, and do not give rise to changes of sign of the IEC. On the 



experimental point of view, this e ffect was confirmed by Bloemen et al. [101 in 



Co/Cu/Co(001) and by Back et al. j|102| in Fe/Cu/Co(001). It has also been confirmed 



theoretically by Nor dstrom et al. 103 , Lang et al. |104| , Drchal et al. [ |105| , and by 
Lee and Chang pi. 



6. Effect of over layer thickness 

A more (at first sight) surprising behavior is the dependence of the IEC on the 
thickness of an external overlayer. From a naive point of view, one might think 
that layers external to the basic ferromagnet/spacer/ferromagnet sandwich should not 
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Table 1. Comparison between the theoretical predictions of Ref. |33[ and 
experimental observations for the oscillation periods of interlayer exchange 
coupling versus spacer thickness. 



spacer 


theory system experiment Ref. 


Cu(lll) 


A = 4.5 AL 


Co/Cu/Co(lll) A ps 5 AL S 


Co/Cu/Co(lll) A ps 6 AL 


92 




Co/Cu/Co(lll) A ps 4.5 AL 


93 




Fe/Cu/Fe(lll) A ps 6 AL 


94 




LI I \J\J ±. J 


Ai = 2.6 AL 
A 2 = 5.9 AL 


Co/Cu/Co(001) A ps 6 AL 


95 




Co/Cu/Co(001) Ai ps 2.6 AL 
A 2 ps 8 AL 


85 




Co/Cu/Co(001) Ai ps 2.7 AL |87| 
A 2 ps 6.1 AL 


Co/Cu/Co(001) Ai re 2.7 AL (6§ 
A 2 ps 5.6 AL 


Fe/Cu/Fe(001) A ps 7.5 AL (74 


Ag(OOl) 


Ai = 2.4 AL 
A 2 = 5.6 AL 


Fe/Ag/Fc(001) A x ps 2.4 AL §§ 
A 2 ps 5.6 AL 


Au(OOl) 


Ai = 2.5 AL 
A 2 = 8.6 AL 


Fe/Au/Fe(001) Ai ps 2 AL 89 
A 2 ps 7 - 8 AL 




Fe/Au/Fe(001) Ai ps 2.5 AL 90, gl) 
A 2 « 8.6 AL 



Table 2. Comparison between the theoretical predictions of Ref. |107J and 
experimental observations for the oscillation periods of interlayer exchange 
coupling versus overlayer thickness. 



overlayer theory 



system 



experiment Ref. 



Cu(001) 


Ai = 2.6 AL 
A 2 = 5.9 AL 


Cu/Co/Cu/Co/Cu(001) A ps 5 AL [ 


108| 


Au(OOl) 


Ai = 2.5 AL 
A 2 = 8.6 AL 


Au/Fe/Au/Fe/Au(001) Ai ps 2.6 AL [ 

A 2 ps 8.0 AL 


109| 



Au(lll) A = 4.8 AL Au/Co/Au/Co/Au(lll) A ps 5 AL |110 | 



influence the interaction between the two ferromagnetic layers. This view is incorrect, 
in particular when the system is covered by an ultrathin protective overlayer. In this 
case, the electrons are able to reach the vacuum barrier, which is a perfectly reflecting 
one, so that strong confinement and interference effects take place in the overlayer, 
which leads to a weak but significant oscillatory variation of the IEC as a function of 
the overlayer thickness [107|. 

This effect, which follows directly from the quantum interference (or quantum size 
effect) mechanism, has been proposed and experimentally confirmed independently 
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by de Vries et al. fllO 



for the Co/Cu/Co(001) system with a Cu(OOl) overlayer, by 
Okuno and Inomata 1 109 for the Fe/Au/Fe(001) system with a Au(OOl) overlayer, 
and by Bounouh et al. |110| for the Co/Au/Co(0001) with a Au(lll) overlayer. In 
both cases, the observed period(s) for the oscillations versus overlayer thickness were 
found to be in good agreement with the theoretically predicted ones. This effect has 
also been confirmed by means of first-principles calculations for the Co/Cu/Co(001) 
system with various types of overlayers [111, 112, 113|. The comparison between 



the periods of oscillations versus overlayer thickness predicted theoretically and those 
observed experimenta lly i s giv e n in Table ^. A more detailed discussion of this effect 



can be found in Ref. [107, 111, 113] 



7. Strength and phase of interlayer exchange coupling 

In contrast with the excellent agreement between theory and experiment which 
is obtained for oscillation periods, the situation for the amplitude and phase of 
oscillations is less satisfactory. According to the theory exposed above, the coupling 
takes the following form in the limit of large spacer thickness (asymptotic limit): 

J = Hj£ sm(q a D + <t> a ). (29) 

a 

Since the coupling constant J has the dimension of an energy per unit area, 
the parameters A a characterizing the coupling strength of the various oscillatory 
components have the dimension of an energy. By taking typical values of the Fermi 
wavevector and velocity, it is easy to see from eq. (^?]) that they are typically of the 
order of 1 to 10 meV. 

Table || presents a comparison of theoretical and experimental values of the 
oscillation amplitude strengths A a , for various systems. || We observe that there 
is a rather strong discrepancy between theory and experiment, but also among 
various theoretical studies. Although the agreement seems to be rather good for the 
Co/Cu(lll)/Co system, more experimental and theoretical data would be required in 
order to know whether the apparent agreement is conclusive or accidental. 



7.1. Co/Cu(001)/Co 

The Co/Cu(001)/Co system is the one which has been most investigated theoretically 
and is considered to be a model system to test the predictions of theory. The 
theoretical results reported in Table |^ correspond to semi-infinite magnetic layers, 
whereas the experimental data have been obtained for magnetic layers of finite 
thickness. As discussed in Section [s] the strength of the coupling varies with 
magnetic layer thickness, which can be a source of discrepancy between theoretical an 
experimental results. An other possible source of discrepancy arise from unavoidable 
imperfections (roughness, intermixing) of the experimental samples. 

Let us first address the short-period oscillatory component (labeled with the 
subscript 1). As discussed in Section ^| above, this component arise from 4 equivalent 
in-plane wavevectors kp located on the T — X high-symmetry line of the two- 
dimensional Brillouin zone |}6| . Since the majority-spin band structure of fee Co 
matches well the one of Cu, one has \r\\ w 0. On the other hand, for minority-spin 

§ Note that, in order to be able to compare various theoretical results with each other, we included 
in the present discussion only the calculations pertaining to semi-infinite magnetic layers. 
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fee Co, there is a local gap in the band structure of symmetry compatible with the 



Cu states 

m 



which leads to total reflection, i.e., 



1. Thus, one has Ar x | m 0.5 



[114, |42| and | Ar x is (almost) independently of Co thickness |105| . The various 
theoretical values for the amplitude A\ listed in Table || agree rather well with each 
other, except the one of Re f. Jll5|| w hich is almost a factor 2 larger than the values 
obtained by other authors [114, 105, 116|. This discrepancy may be due to an error 
in the estimation of the curvature radius Hi of the Fermi surface, and of the Fermi 
velocity vj_i, which are quite tricky to obtain accurately for k|n. 

Turning now to the comparison between theory and experiment, we notice that 
the calculated values of A\ are considerably larger than the measured ones. At least 
two reasons can contribute to this discrepancy. The first one is the effect of interface 
roughness, which generally tends to reduce the amplitude of the coupling oscillations 
[j33| ; this effect is particularly pronounced for short-period oscillatory components, as is 
indeed confirmed experimentally |S7|| . The second reason is of intrinsic character: the 
theoretical values of A% given in Table || correspond to the asymptotic limit, whereas 
the experimental data have been obt aine d for spacer thicknesses below 15 A L. A s 
appears clearly from Fig. 6a of Ref. 105 and from Fig. 13 (bottom) of Ref. [116 1, 
the asymptotic regime is attained only for thicknesses above 20 to 40 AL; below, the 
envelop of the oscillations deviates significantly from a D~ 2 behavior, and the apparent 
amplitude in the range relevant to experiments is typically a factor 2 smaller than the 
asymptotic amplitude. This preasymptotic correction is attributed to a strong energy 
dependence of r[ [116]. 

Let us now discuss the long-period oscillatory component. As appears from Table 
||, the situation is quite confusing: no t only th e various theoretical results disagree 
with each other, but some of them 
as compared to the experiment 




15, 1 16 1 underestimate the coupling strength 
a fact which cannot be explained by the 



^0] |S5| 

effect of roughness or interdiffusion. 

The long-period oscillatory component arises from the center T of the two- 



dimensional Brillouin zone. Here again, for the same reason as above, one has \r 2 \ ~ 0. 
The minority-spin reflection coefficient, on the other hand, is considerably smaller 
than for the short-period oscillation, and one has 
M 



0.15 pa], so that | Ar 2 | « 0.05 



114 . 



This very small spin-dependent confinement explains the very small values 
of A2 obtained by the authors which rely on the asympt otic expression (^?j) obtained 
from the stationary phase approximation 114 115 , [ll6| . However, as seen from Fig. 
2 of Ref. [118] and Fig. 2 of Ref. |115|, r\ increases very strongly with k|| and full 
reflection is reached at a distance 0.1 x n/a from T; indeed, the low reflectivity arises 
only in a narrow window around T. As discussed in Ref. [ 1 19 1 , this gives rise to a 
strong preasymptotic correction, and explains why the stationary phase approximation 
approximation yields an underestimated value of A2- On the other hand, if the kii 
integration is performed without using the stationary phase approximation, as was 
done in Ref. |105| , a much higher value of A2 is obtained; the latter is larger than 
the experimental one 



p6| by a factor 2.5, which seems plausible in view of the 
effect of roughness and interdiffusion. 

Our knowledge of the phase of the oscillations is much more restricted as this 
aspect of the problem has attracted little attention so far, with the notable exception 
of the work of Weber et al. f87fl . On general grounds, in the case of total reflection (as 
is the case for rij), one expects the phase to vary with magnetic layer thickness and/or 
with chemical nature of the magnetic layer; conversely, for a case weak confinement 
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Table 3. Comparison between the theoretical predictions and experimental 
observations for the oscillation amplitudes A a of interlayer exchange coupling 
versus spacer thickness; for Cu(OOl) and Au(OOl) spacers, A\ and A2 correspond, 
respectively, to the short-period and long-period oscillations. 



system 



theory 



Ref. experiment 



Ref. 



Co/Cu(lll)/Co 



.4 
A 



3.7 
4.2 



mcV 
mcV 



[114| 
[115| 



.4 
.4 
.4 



7.6 meV 

3.4 mcV 

2.5 meV 



117 



Co/Cu(001)/Co 





a 42 meV 


A 2 F 


a 0.13 meV 


A x f 


a 72 meV 


A 2 f 


a 0.75 meV 


A X F 


a 35 meV 


A 2 F 


a 3.5 meV 


A X F 


a 35 meV 


A 2 F 


a 0.035 mcV 


A X F 


a 12.5 meV 


A 2 F 


a 6.9 mcV 



111 



1 115| 



[1051 



A 2 



1.6 meV 
1.4 meV 



Fe/Au(001)/Fe 



(1151 



Ai « 8.1 meV 
A 2 sa 1.1 meV 



(as is the case for rJ;), one expects the phase to be almost invariant ][36|. These general 
trends are indeed confirmed experimentally by Weber et al. |$7J . 

7.2. Fc/Au(001/Fe 

The system Fe/Au(001)/Fe is actually an excellent system for a quantitative test of 
the theory. This is due to the excellent lattice matching between Au and bec Fe 
(with rotating the cubic axes of the latter by 45°), and to the availability of extremely 
smooth Fe substrates (whiskers) |9C], |9l| . 

In contrast to the Co/Cu(001) case discussed above, for Fe/Au(001), one has 
total reflection of minority-spin electrons both at (short-period oscillation) and 
k|| 2 (long-period oscillation), and \r+\ is almost independent of k | around these points, 
as appears clearly from Fig. 1 of Ref. ][l2|. Therefore, the associated preasymptotic 
correction should not be very strong. 

Indeed, as is seen from Table pi the predicted amplitudes are quite large, both 



for the short-period and long-period oscillatory components [115]. These predictions 
are fairly well confirmed by state-of-the-art experimental studies |9l[], although the 
predicted amplitude of the long-period component is too large by a factor 6. 

Clearly, even for this almost ideal system, further work is required to achieve a 
satisfactory quantitative agreement between theory and experiment. 

8. Concluding remarks 

As has been discussed in detail in this review, there is a great deal of experimental 
evidence that the mechanism of quantum confinement presented above is actually 
the appropriate one to explain the phenomenon of oscillatory interlayer exchange 
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coupling. This mechanism is entirely based upon a picture of independent electrons. 
This may seem paradoxal at first sight, in view of the fact that exchange interactions 
are ultimately due to the Coulomb interaction between electrons. In fact, this 
independent-electron picture can be be justified theoretically and is based upon the 
"magnetic force theorem." A thorough discussion of this fundamental (but somehow 
technical) aspect of the problem is given elsewhere [120, 12 1 1 . 

In spite of the successes encountered by the quantum confinement mechanism, a 
number of questions remain to be clarify for a full understanding of the phenomenon. 
In particular, one needs to assess in a more quantitative manner than has been done so 
far the validity of the asy mpto tic expression ((2^); a first attempt towards addressing 
this issue is given in Ref. [ 1 1 9 1 . 
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